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ABSTRACT 


Avdasieticstoneot the problem of the irreducibility of 
polynomials in the ring of integral polynomials establishes 
the framework of the research. A transformational scheme 
is postulated to facilitate investigation of the problem. 
The coherency of the scheme is detailed and the necessary 
Computational techniques developed. 

To determine the efficacy of the transformational scheme, 
tie Specification and collection of appropriate sets of data 
are discussed. The transformational scheme is then applied 
to the data aif UicmkesultS are Cabwlated and discussed. 

The conclusion is drawn that a transformational scheme 
memceucerw) tool fer the investigation of the irreducibility 
of this type of polynomial and some hypotheses as to the 


forrnement and extension of the technique are stated. 
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I. INTRODUCTION 


fepolynenaal is said to be reducible in the ring of 
Dot yiemradlice wWisen integer coefficients, called Z[x], if it 
is expressible as the product of polynomials in Z[x]. A 
polynemaal not reducible as defined above is irreducible in 
Zales. 

The polynomials under consideration in this research 
were those polynomials in Z[x] and reducibility over the 
Pie worteimcem@ers, Z. The general formeof this type of poly- 


nombal will be stated as follows: 
p(x) = a,x + a,_4X to see + aX + ap (1) 


investigation into the irreducibality aspects of 
polynomials leads to two primary avenues of research. The 
Hivest MaiacOr May eatbenpt to detenumame whether a given poly- 
nomial is irreducible or not. The determination is made as 
ume reswit of some analysis of the polynomial. The analysis 
Mey weconscise of identifying some pattern among the coefficients 
an some predetermined relatiemship among the coefficients, or 
some other type of analysis. | 

Tieeprenltem tien becomes one of finding the appropriate 
Speeciiiedqeron for the given polynomial or class of polynomials 
which will resolve the irreducibility of the polynomial or 
elluge Smee monials Under Imvestigation. These specifica- 
tions for individual polynomials or classes of polynomials 


imme as “mereducibility criteria’. Several criteria are 





known and apply to many different polynomials (Dorwart, 
ies) However, no gemeral arreducibility criterion has 
been found which roll apply to all polynomials in Z[x]. 

Alternatively, the investigator may attempt to extract 
the factors of a given polynomial directly without first 
deciding the question of irreducibility. If the process 
results in the extraction of factors then the polynomial is 
reducible; otherwise, it is irreducible. This approach has 
resulted in several algorithms which extract certain types 
Gt taetorms (Knuth, 1969). The Kronecker method (Dorwart, 
WSS) aise a general solution to the problem of factor ex- 
tmactrom of polynomials in Z[x] but it is awkward to apply 
and the calculations involved are prohibitive in length so 
Guat the method iseof little uté#hity for practt#eal work. 

An illustration of this point will be discussed ihawer. 

Thsss, both paths of attack into the area of irreducible 
polynomials in Z[x] have been less than completely success- 
ful. A general irreducibility critericn, if one exists, has 
not been identified and the general solution to factor ex- | 
tmaaetion is pPeohibitive in length» even with the speed of 


Calculation provided by the electronic computer. 





Il. TRANSFORMATIONAL SCHEME 


Paeeer LTeviewine Several Studies in the area of irre- 
Pucrple polynomials in Zixj, 1t was decided that rather than 
attempt to solve the problem of irreducibility directly for 
a given polynomial, perhaps it would be more fruitful to 
attempt to facilitate the solution by somehow transforming 
the given polynomial into a polynomial for which the irre- 
ducibility had previously been decided, or would be easier 
to decide. 

This decision as to a method of attack led to a search 
for transformations which could be used to, in some sense, 
simplify the polynomials under investigation. To be of any 
aeteerOrenais work, the transformations #rould haere ve be 
‘such that the irreducibility qualities of a polynomial would 
be preserved across the transformation. 

iWewt1rst transtommation which suggested itself was to 
Beaistorn am arbitrary polynomial with integer coefficients, 
say f(x), to a primitive polynomial, i.e., a polynomial, say 
g(x), such that the greatest common divisor of the coeffi- 
Gmemts 1S Unity. The irreducibility of f(x) is preserved 
Srmecesii%) 1S the product of g(x) and some integer k, which 
iommme gcd of the coefficients of f(x). Clearly, if £(%) 
is irreducible, then so is g(x) and vice versa. 

Given @ primitive polynomial, it may be further sam- 
plified by transforming it to a polynomial whose leading 
coefficient is unity, known as a monic polynomial. Consider 


a primitive polynomial with integer coefficients 





xX Te © @ eee Fe xs a 


i 1 ge, 


where a #1. Then a transformation [3] called M, such that 


M 
-] 
p(x) —> at o( | = q(x) (3) 
n 
and 
Bo) 
= x" n-1 ee n-2 n-1 
ques): x + a. 4X + a,4a,_2% + aes a,x + a a 


produces a monic polynomial with integer coefficients (MPWIC). 
fire question arises as to the preservation of irreducibility 
across this transformation. 

Assume that p(x) is irreducible in Z[x] and suppose that 
q(x) is reducible in Z[x]. Note that q(x) is monic, hence 


peamrcive, then 











Giese) = ac (X)> 19) (5) 
ae if x). exo OX) (6) 
[E+ * k(a,x) r(a,x). (7) 
Nn a 
n 
substituting a,x Ol eX: 

£(x) + cr ) ay r(a_x) (8) 

a 

nN 


Memecn tix) = k'(x) r'(x) where k'(x) and r'(x) are in Q{[x]. 
RwteEsGaiss'’ theorem states that £(x) 1s reducible in Q[x] if 


icmonmly Tf 1t 1s reducible in Zix]> ~Yinerefore we have a 





Eomereaatection and must comelude that if p(x) is irreducible 
then q(x) is aliso irreducible. This result establishes the 
feescrvation Of irreducibility and gives the result that if 


mee 1S a primitive irreducible polynomial with integer co- 


n-l 


efficients then q(x) = an 


p(x/a,) is a monic irreducible 
peenomtal with integer coefficients (MIPWIC). Since both 
the reduction to primitivity and the transformation to mon- 
icity preserve irreducibility, it follows that any polyno- 
‘meat may be transformed to a MPWIC with irreducibility 
properties identical to the original polynomial. Thus the 
only polynomials to be wimmestigated fon arreducibpility sane 
monic memeaeni as With integer coefficients. 

Attempts were made to discover other transformations 
which would further Simplify the MPWiCS's, while preserving 
Meontcity, irreducibility, and eet GOcurrercmes.. Meme 


transformation considered was 


f(x) ——— £(xtk) : most (9) 

where if 
wx) = x + ea t+ see + ax + ap (10) 
then (11) 


£(xtk) = (xtk)™ + a,j (xtk)™ D+ eee + aj (xtk) + ap. 
Clearly, the transformation preserves monicity and integral 
Geemticienes., AS to the preservation of irreducibility, 
SMmompese that £(x) is irreducible and f(x+k) is reducible. 


omsider t(x+ky = g(x) h(x). Replacing x by x-k we get™ the 


expression 





f(x) 


ll 


~e(x-k) h(x-k) (12) 


£(x) = gi(x) h'(x) (13) 
wemeress (x) and h'({x) are in’Z{/x}. Therefore £(x+k) = 
pues) nix) implies that f(x) = ge{x-k) h(x-k) which is a con- 
tradiction of the assumption that f(x) Sse cues 
Hemec, irreducibility 1s preserved across the transformation. 
This transformation will be subsequently referred to as the 


'T' transformation. Other notation used to describe the 


transtormation will be 


T 
E(x) —“) g(x) and fT, = g. (14) 


The next transformation considered will be known as the 
'D' transformation in that its effect on a polynomial is a 


@ilation of the axes. The D transformation is specified by 


Go) =) ie ele , n= deg (f(x)). GUS) 


Using similar notation as for the T transformation, other 


forms are 


D 
ec eS g(x) and £D, = g. 


fie crrect Of this transformation is given by the following 


peoomess@an. If 


tees tl ipo -_. 
Ee). =ex. + a4 + + a,X + ap (16) 
then (ie) 
Sic £( + eee Sa i. eaee Ot ce a : ze 
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Clearly, monicity is preserved. If the number b is an in- 
teger then the transformation preserves integral coeffi- 
eeemes. However, this number may be the reciprocal of an 
integer if the coefficients of f(x) are such that b" divides 
2a for each coefficient of f(x). The preservation of ir- 
reducibility is proven by methods analogous to the monicity 
transformation using Gauss' theorem. 

A Ghieed tramsitormation, called the X transformation, 
was considered. This transformation consisted of replacing 
the variable in f(x) by its reciprocal and multiplying the 
resulting polynomial by x", This transformation was shown 
to preserve aaeeual SOen herent cmda te br Palienakioenion but 
preserved monicity only in the case where the last coeffi- 
ememt OF 1{X) waS UNIByY since wes cifect on a polynomial is 
aaevciscwine Order of the coefficients. Thus the X trans- 
formation was eliminated from further consideration as being 
too limited in application for the purposes of this work. 

Vemmmeethiatecach different valwe of k or b for the T 
ange D transtormations will produce a unique irreducible 
polynomial it is clear that a single irreducible polynomial 
Generates an infinite number of ee Cree Orne cures 
This leads to the idea of ‘equivalence classes' of polyno- 
mials with an equivalence class being defined as a set of 
irreducible polynomials of like degree, each of which could 
be produced from any other member of the set through a 
properly chosen sequence of D or T transformations. 

several relationships were developed which relate the 


effects of sequential series of the two transformations. 


iil 





T T 
f(x) —K4 £(x+k1) —K2y £(x+k1+k2) (18) 


ie, ot kK' = ki + ko, then 


Ty 
egw F(xek') = E(x¥ky+k,). (19) 


Thus, the effect of any pair of successive T transformations 
can be produced by a single T transformation with a properly 
chosen value of k. Through veneated Sade ait Onsi Od Ea@tia - 
tion (19), it is clear that any finite sequence of T trans- 


formations can be replaced by a single T transformation. 








Simplarly , 
D D 
£(x) ae Pe ees mes be be ee (20) 
1 *|d, 1 3 5 
but, if b, = b,b,, then 
Bes X _ Aaa oe 
f(x) ——5 bo £ a = b? be {5}. 21) 





iimerefore, the effect of a@ny two sticcessive D transformations 
can be achieved by a single D transformation through the | 
Peeper GCnolce of the tranStormation factor. 

Now, consider successive transformations of different 
ipes, iee., a DD transformation followed by a T transforma- 
tion, or vice versa. Consider the TL Dy pair in the following 
expression : 

qT 7 a n x 
f(x) — > f£(x+k) —> b ‘(i | (22) 


But, after rationalization, the result is 


NGA 








3 


D T 
b n x k n 
f(x) ——> wv” £(t]) ——> dv” 2% (24) 


In this case, simple algebraic manipulation yields the fol- 
lowing expression: 


' Dix. cn fx . k 
(x) —2S vo [+ El, (25) 


which is the same result as would be affected by the Te pp 


u@aiistormation pair, From the establishment of the two pair 


relationships 
ee Dale (26) 
and 
Dat, = ie ; (27) 
a 


alone with Equations (19) and (21), it follows that any 
finite sequence of D and T transformations may be collapsed 
into a DT or TD pair of transformations, by first collapsing 
the runs of like transformations and then reversing unlike 
pairs to produce successive like pairs which may again be 
collapsed into single transformations. The reduction pro- 
cess may be continued until the effect of the sequence is 
represented by, at most, a TD or DT pair of transformations. 
Sumlede transtommations, either T or D, are included in 
Gee case GOrwpairs since, from the definition of the trans- 
formations, it follows that a To ane D, transformation returns 
the original polynomial. Thus, any single transformation may 


be considered to be part of a DT or TD pair, where the other 


ite 





PaaisieGiiiattonm im the pair is the identity transformation 
for the unspecified type of transformation. Note also that 
each type of transformation has its inverse. The inverse 
Ob chic TL transformation is T_, and the inverse of the Dd. 
transformation is D : 
iy ai 
Summarizing, a transformational scheme has been devel- 
oped, and an algebra of the operators established, which 


will ‘permit the manipulation of the polynomials under in- 


vestigation in this research. 
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Ili. DEVELOPMENT OF ALGORITHMS FOR TRANSFORMATIONAL 
OPERATIONS 

The ne oe on ad scheme described in Section II 
establishes a coherent framework for manipulation of poly- 
nomials. However, in order to actually perform these opera- 
tions on a digital computer, a series of algorithms must be 
developed. Further, to be of more than nominal value, the 
pueorrthnms have to be general, i.e., perform the desired 
operation on a monic polynomial with integer coefficients 
of degree n where n is a positive integer greater than or 
Gamal to 2. 

The first operation to be considered is the single D 
transformation. By defamiita on 


D 





ce 
£(x) ——e b” £(=| . (28) 
Let 
= n n-1 ee 6 
f(x) = a,x + a,_4x + + a,X + a) (29) 
then 
(30) 
Boece = ee nen + q ce, perc | aa 
b : n n igs: ipl) 1\b 0 
b b 
7 n Tite ek Tae igh= n 
=a x + n-4>* 1: 25 a,b x ab ($2) 
n 
= y Wee, | oe ie (32) 
1=0 


Thus, Equation (32) establishes an algorithm for computing 


the coefficients of b” £(x/b)seamen the coefficients of 1(%), 


1S 





the degree n, and the value of b. Programming of the al- 
gorithm is straightforward and will not be discussed fur- 
Emer . 

An algorithm for computing the coefficients of the re- 
sult of a T transformation was more difficult since, except 
for the leading coefficient, there is no simple relationship 
between the coefficients of the original polynomial and the 
transformed polynomial. The definition of the T transfor- 


itierOn gives the expression 


ik 
k 
f(x) —— f£f(xtk) . (53) 
et 
f(x) =a ae OT lil + + a,x + ap (34) 
then 
5) 


f(xt+k) = a (x+k)" th OSS) Sie enone a, (xt+k) + ap. 


For an algorithmic solution, consider Taylor's formula, as 


Piven py the expression 





; (n) on 
ees 2) ee), Pe) ra) 56) 


Pugetreuting xtk for x and letting k=a, Taylor's formula 


becomes 


(n) n 

aes AK), Eee, ge (ax! (37) 
| 0! 1! n! 

Vatetewrelas tae Interesting result that the coerficient, bs 


of f(x+k) is given by 
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1 1! (38) 
Since 
f(k) = ak" 1 hal 1: ine ice + a,k + ap (39) 
then 
(40) 
es ee tae (n= 1) ka tae (n-=2)k" 24 >>. aya 
n n-l n-2 1 
£6™°™ (x) = a (n)(m-1) +++ (n-141)K™ + a (n-1) (n-2) 
41) 
m-1 ( 
(n-1)k ie eC ee See 
(42) 


e-VWeky = a (n)(n-1) +++ (2) k + ay (m-1)(n-2) ++ (CD) 


(0) (k) = a(n) (itl) Tere (1). (239) 


Statins the general term-of the sequence in series notation 
piwes the expression 
Mm ° 
é m-1 
Cien) . (erin) 
: (x) 4an-i  (m-ayt 2 
i=0 


Substituting Equation (44) into Equation (39) for ft ae 


J i 
n oe . =: 
| : Coa elles 
i 


n-i Gane 
£(x+k) = ee 
mi] )! 


a0 





Equation (45) provides an algorithm for computing the coef- 
ficients of f£(x+k) given the coefficients of f(x), the 

degree n- of cowed a value of k. Again, as in the D 
transformation algorithm, the programming is straightforward 
except for the computation of the necessary factorials, which 
Can be handled by an iterative or recursive function. 

The development of the algorithm for the T and D trans- 
formations permits the transformation by computer of any 
polynomial by any value of k and any non-zero value of a. 
However, since the primary g@al of the research, hence, 
transformational scheme, is to reduce sets of polynomials 
into classes rather than genemase classes of polynomials, 
investigation of methods of determining whether a transfor- 
mation exists between two arbitrary polynomials of like 
demmee 15 of greater interest. | 

The cases of two polynomials differing by a single trans- 
formation or by a sequence of like transformations will not 
be discussed separately as the relationship of a single trans- 
formation to a pair of transformations was discussed in see: 
tion II and an algorithm for@a pair will inciude the cases 
of single transformations. As has previously been establish- 
ed, any sequence of T and D transformations can readily be 
Gomlansed antema DI or TD pair and, further, a TD pair can 
be replaced by a DT pair of transformations and vice versa. 
Thus, any algorithm to determine the transformational rela- 
tionship between two polynomials of like degree need only 


determine if they are related by, at most, a pair of unlike 


18 





tetoronnarions, Rather than test for both types of pairs 
an algorithm which determines only a DT or TD relationship 
Pemsuti@etrent since the existence of one pair implies the 
existence of the other pair between the same two polyno- 
mials. 

On examination of Equations (26) and (27), which es- 
tablished the relationship of the two pairs for any two 
Daimemals, motesthet mm the DI pair the total translation 
mrecor 25 Cither k or ka. Hence, if k and a are both in- 
@epers, the total translational factor of the pair is also 
mmeitercr, Also, from the definition of the D transforma- 
Een, Mete that the factor, a, must be either an integer or 
eiemrceiprocal Of an integer if the coefficients of the 
Paroroermea polynomial are to be intégers. From the state- 
ment of the inverse of the D transformation, it is to be 
noted that if the om! transformation pair, which- 
ever exists, from f(x) to g(x) is such that a is not an in- 
teger, then the DoT, or De pair, respectively, from g(x) 
Peeieoyewniie be such that a is an integer. Then, since k 
must always be an integer, between any two polynomials which 
are related through a sequence of transformations, there 
eeeeotcmie)) pair such that the T factor and the D factor 


are both integers. That is, of the four relationships 


Dat 
21 Ii D Ty 


oe ———— fC) Cae E(x) ee g(x) (46) 


Pe ka 


Cl Seen. 6X) Gee f(x) ———> g(x) (47) 


Is 





1D eel i TD 


aval nrg 
eo) a ieee (x) ——— ££) (48) 
la Doty 


Oa 7 eee (x) = ——> f(x) (49) 


restated here for clarity, at least one must be such that 
both transformational factors are integers in the DT pair. 
Hence, the algorithm has only to determine whether or 
ferme at most a DI pair with integral factors exists between 
twe polynomials to completely decide whether they are rela- 


ted by any sequence of transformations. 





Let 
_ n n-l . 
f(x) = bx + b-1% + w byx - bo (50) 
and 
x oT n-l : 
g(x) = c,x + CL_4Xx > * C)X * Cy ese) 
be two polynomials such that 
vavk n Se) S 
f(x) —-— a f = 7 Palle) an. (52) 
Subst tution in Equation (52) gives the expression 
n Nace Tl xtk n-1 xt ke = 
aa ee a a veer) eon 
n ne . 
c\x + ¢,_4x ‘ + C)X + Cy (53) 


By repeated differentiation of both sides, we obtain 


Hee) +k =, fn at 1 
cg a = Dy | Z hl 

















— 


m + bi.z(n-1)! 
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di x+k = Te 2. ; 
“n-2 e | + bi-2(n-2)! mes 24) oo C1 {a ee ox 
i C,-2(n-2)! = g ee Coe (49) 
and 
Meett- 1) ei xtk] _ nn 1 xt r il a 
a is ae bn! [4] a | ee ha 1" 





c,(n!) oe c,.z,-1)! - gM™-D) xy CS 


mmee Fauations (54) and (55) are true for all values of x, 
they are certainly true for x = 0. Simplifying and evaluat- 


ing at x = 0, it follows that 


2 
2 Did Wik 2 k _ 
. », | ns (&| ia bp, -pea-w hs B,-2tnr2yt = 
c_4(n-2)! | (56) 
and 
as n'( 4) 2 als. (asin e_ . Gresley (57) 
n a n-l ° n-l ite , 


Hurther simplification produces a system of two equations in 


two unknowns, k and a. 


(58) 
:, . ? 
kK? [ban én-3)| oils [2b,-180-0)] 3 | 2.28 = 2-2] = 0 
and 
b kn bere 2 = | Coery foo) 


Somme Equation (59) for a in terms of k and substituting 


for a in Equation (58), an equation in k is produced which 


rag 





Memecsaly solvable by the use of the quadratic formula: 


2 2 2? 
k | b,mbp 1 Cl-m + 2bon ba. | 


2 2 
+ k | emmepal(ba = Zbaebs - 2y-2P aot] (60) 


4 | a 2b ena | o 

PaoesoOlulrTonmo: Equation (60) is a non-negative in- 
teger, Equation (59) is solved with this value of k. If a 
ise also an integer, then the DoT, EransStoOrmatlon 1S per 
hermea om £{X) Wsing the algorithms previously developed 
Siemene result is compared to g(x). If £(x) DT, = g(x) 
tiem ehRere iS asrélationship between f(x) and g(x) and the 
factors of the transformation are k and a from Equations 
(59) and (60). If neither solution of Equation (60) is 
Satisfactory then Equation (60) is modified by substituting 
bs Low Cs and C5 One bs and ene result is solved in the same 
manner. This interchange technique allows testing for an 
iateeeral Palr Of transtormations from g(x) to f(x). If 
g(x)D,T, = £(x) for some values k and a then again there is 
ala eLoneanimmamednat 1s.specified. If not, there exists 
no transformational relationship between £(x) and g(x). 

The implementation of the algorithm is straightforward 
except for the special cases of b,j and/or c,_, being zero 


Since Equation (59) cannot be solved in these cases. The 


“ZL 





Beenmnique utilized in these special cases is to solve Equa- 
tion (59) for k instead of a and then substitute this solu- 
tion into Equation (60) giving the following system of 


equations: 


2 
a [by Cy." 2)(@-2) + 2b» 2? | 


(61) 
a | n-1 €n-12) C2) 23 | eA Fituneeds) = 2b aon 2 | = 3 


anid 


b_nk + ab = Cc 
n n 


i (62) 


n-l 


Eamaetem (62) 1S always solvable since the degree nm 15 
femal to or greater than 2 and De = 1 due to the monicity 


eemstraint. Clearly, if b = 0 then k = 0 but this 


n-l 7 ene 


gs expected since, from Equation (38), 


Z ¢ (4) (k) 
eat Le) 
hence, 
nt n! ve 
Sel = bo a k + og : : (64) 








Tits, since DD — lente Condition kK = QO 1s the only one 
which, along with the prior condition that a,_, = 0, will 


allow b “m0, Imerefore, if both an] and b__, are equal 


tue 
to zero, the transformations, if any, between f(x) and g(x) 
must be solely of the D type since k = Q. This possibility 


is tested by solving for a D transformation factor from bo 


ZS 





and Cg: Here 1S, af g(x) = a £(x/a) ener Cy = a''by- Solv- 
ie £Or integral values of a and then calculating the trial 
polymomial as in the general case allows the testing of 
Gers Condition. 

With the development of the algorithm for testing f(x) 
and g(x) pairs, the necessary algorithms have been developed 
Eersperforming cither the D or T transformations on @ given 


polynomial or determining the transformational relationship 


between two polynomials. 
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Iv soeeClPICATION AND COLLECTION OF DATA BASE 


aortvccmmenesutility amdeetticacy of the developed trans- 
formational scheme and the associated algorithms as a tech- 
image for investigation of irreducible polynomials in Z[x], 

a representative data base of irreducible polynomials was 
necessary. To provide a large but finite number of poly- 
nomials from which to select the irreducible polynomials, 
certain conditions were specified and all polynomials meet- 
ing these conditions were selected and then screened for ir- 
meaucibility. The conditions for selection were as follows: 

(a) £(x) is a monic polynomial with integer 

@oerticients 

ieee etx) 15 an Z2 ©ror kK in Z 

(c) the absolute value of f(k) is less than or 

equal to 5° for k = 0,1,...;,;n where n is the 
degree of f(x). 

Gemaimions (b) and {c) imply a grid of points in the x-y 
Puaemesewameh is represented graphically in Fig’ 1.° The data base 
zs the set of all polynomials satisfying the three condi- 
Puamcmspecitied above. Note that monicity, by definition, 
Pequiregs that the léading coefficient, an be equal to l, 
rather than |a,| = 1. But if a polynomial, say f(x), exists 


Pileeirnateie saticeuares all the conditions for inclusion ex- 


Cee thiat aa -1, then there exists another polynomial say 
s(x), such that f(x) = -g(x) and g(x) will be included in 
Gg@emaata base. That is, eacheof the monic polynomials in- 


aicdecaein the data set will have a reflection across the 


Z5 





£(b.) 





Penge | . 


x-axis which lies on the degree grid but is not included in 
mHedeka base. But the factorization or irreducibility of 
the reflected polynomial is iden Glee, thigteo: ts re— 
Pieeeionewith the additional factor of -1. Hence, the ir- 
reducibility or factorization of both a polynomial, f(x), 
‘and its reflection, -f(x), are represented by the inclusion 
of f(x) in the data base. | 

The actual calculation of all the polynomials for a 
given degree grid can be accomplished by use of the LaGrange 


interpolation formula 


Oe 
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where bs are the abcissae and the f(b.) are the ordinates 
of the n+l points for the summation. After the calculation, 
the polynomial can be screened for monicity and integral 
coefficients. Clearly, this technique requires that a 
LaGrangian polynomial be calculated for every possible set 
of f(b.) over the grid. Note that the number of sets of 


mimetional values has the order of yyntl 


Where ne 1s ‘the 
maximum degree of the LaGrangian polynomials. Hence, the 
number of polynomials to be calculated becomes extremely 
marge for grids vot modemate degree. The calculations for 
each polynomial consist of the actual calculation of the 
polynomial and en testing the polynomial for monicity and 
antegral coefficients. After programming the LaGrange for- 
mula and executing test runs with Limited grids, 1t was de- 
termined that approximately 780 microseconds were required 
momeeemdch fifth degree polynomial. Extrapolating to the full 
pemeetor the fifth degree, this technique would have re- 
ouerred just under 24 hours of central processing unit time 
Oman 16M 360/67 for the processing of the 1,771,561 possi- 
ble sets y. pommuse Thus an exhaustive search technique, 
which is analagous to the Kronecker method of factorization, 
Msepronibrerve in terms of time required. 

Examination of the grid and LaGrangian formula, along 
with the statement of the problem, offers several areas for 
reduction of the number of polynomials to be calculated. 

The previously mentioned case of reflected pairs allows elin- 


memeron of half the sets of functional values. Thus, if the 


Pee) | 





Premconareron that £(0) 2 0 is imposed then one member of 
each reflected pair is eliminated and the number of poly- 
Wemtals to calculate is halved. Additionally, since the 
@eai 1s to find irreducible polynomials, any set of points 
such that at least one f(b.) = 0 may be eliminated since 
moreany kK, £(k) = 0 implies that k is a root and hence (x-k) 
ts a factor of f(x). Therefore, f(x) is reducible. 
Considering the LaGrange formula and actual calcula- 
Peomieot the polynomials, note that the requirement that the 


pelymeomials be monic imposes constraints on the possible 





values of f(b, ). The LaGrange formula may be written in 
the form 
n 
Lk) bso 
1=0 
where 
n 
Yas) | | (x-j). (67) 
i 
j7k 
Then 
£§) (9) Y 60 a neh (68) 
K=0 "k 


but, since Yy is monic by construction 


yt") (9) = nt | (69) 


hence 
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nN 


fiM) (9) = 4 £ (k) y tky j (70) 
K=0 : 


Eyecefiaition, ¥_(k) = k(k-1)...(1)(-1).-.-(-(n-k)) thus 
Y,(k) = k! ay ee. | (cab) 


Substituting Equation (71) into Equation (70), we have the 


expression 
n 
(n) = n! 
f LOja= ». £(k) a. (72) 
k=0 ee ta) (ike) 
or 
n 


_ aysinignds oe 


(73) 


ray 
© 


But £61) (9) = £{1) (x) which is just the leading coefficient 
Ge £(X). Since the goal is to iewe monic polynomials then 
@m@ly those sets of points which satysfy 

n 


1 = pi £(k){R} (-1) 7% | (74) 


k=0 


can possibly be functional rile of a monic polynomial. 
Thus, only sets of points which satisfy Equation (74) need 
be considered and the LaGrangian polynomials need be calcu- 
Patcumuon sony these sets of points. 

Mmeoraer to illustrate the savings effeeted through the 


use of the condition established by Equation (74), consider 


Zo 





the 1408 minutes of CPU time necessary to calculate and test 
all the sets for the fifth degree grid. Elimination of the 
memlected pairs sets and those sets containing a zero still 
leaves 500,000 sets of points to consider which would re- 
quire approximately 398 minutes of CPU time. The conditions 
of f(b, ) # 0 and Equation (74) were programmed and these 
500,000 sets were tested. In this case 551 sets were found 
mi apprexanately 38 seconds of CPU time. The S52 LaGrange 
Porynonials were then calculated and tested in another 32.5 
peeends fer a tetal of /0.5 seconds compared to the 398 
mmmces required for the exhaustive search technique de- 
scribed above. | 

The same conditions were applied to the grids for 
degree 2,3, and 4 and a limit of +5 for each case. The 


Pesults are summarized below: 


TABLE I 
NUMBER OF 
SETS SATISFYING MPWIC!S | 
LIMIT X DEGREE EQUATION (74) FROM LA- 
Ge TOTAL SETS AND f(b.) #4 0 GRANGE 
_ FORMULA 
SO 500 44 44 
oes 5000 7 256 Wee 
ex 4 50000 1067 Lay 
Bes 500000 551 0 


The monic polynomials calculated for each grid (see 


Appendix A) are a superset of the irreducible MPWIC's for 
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e@en grid since they are the set of all the monic poly- 
nomials with integer coefficients except those which con- 
jaan a linear factor of (x-b;) Since those polynomials have 
already been shown to be reducible. However, polynomials 
men linear factors of the form (x+k) have not been elimin- 
ated. Thus each set of polynomials was further screened by 
a synthetic division program to eliminate those polynomials 
meen Contained this type of linear factor. Since the ab- 
Solute value of the constant term of any linear factor of 

a polynomial must be less than or equal to the absolute 
value of the constant term of the polynomial if the poly- 
fembal is reducible in Z[{x|, the synthetic division program 
tested those values for each polynomial under test. This 
met, combined wath the previows conditions imposed upon 
mae sets of polynomials, constituted a sufficient condition 
fer irreducibility for the MPWIC's ofedegree 2 and 3 since 
if MPWIC's of these degrees are reducible in Z[{x], they 
must contain a linear factor. However, those polynomials 
omecegrees > 4omay not contain a linear factor and still 

Pere dicipwe in Z[x!] sas the product of two imreducibles poly- 
nomials of degree eh 


Suppose that 


ee ee) 1X) (75) 


then for any value k, f(k) is the product of g(x) and h(x), 
ic. , 2k) Ws divisible by g(k) and h(k) for any value of k. 


fae x<jywandsin(x) améwallaneeiels “then for am integral 


By 





Memeo oOr k, £(K), @(K), ama h(k) are in Z. Also, if 

£(k) # 0, then g(k) # 0 and h(k) # 0. Now if an integer is 
me produec of two ormmore integers, none of which is zero, 
then the absolute value of each of the factors of the pro- 
Gm@et 15 fess than or equal to that of the product. That is, 


merry xy = g(x) h(x) then” for any value of x 


|£(x) [>1e¢x) [21 (76) 


and 
|£(x) |>|hGo [21 . (77) 


Applying these results to the polynomials under discussion, 
it follows that if f(x) is a monic polynomial with integer 
coefficients of degree > 4, |£(k)|<5, and f(x) has no lin- 


, then its factors must 


ear factors but is reducible in Z[x] 
be polynomials over a grid of lesser degree with bound equal 
to §. In the case of the monic polynomials, f(x), of degree 
M@ with no linear factors and such that |f(k)[<5 for k = 
we. .,4 then any possible irreducible quadratic factors 
mst be ancitided in the set of Pe auGivac quadratacs Yal- 
ready determined. Thtise. Coal Givwicwon von bthe atarwices wath 
no linear factors by each of the iserechie ple quadratics will 
Poros mene 11reducibility of the quartic in question. 
The actual technique employed was the testing of each co- 
Perieremt. ot the quartic successively during the division 
so that the full division occurred only if the quadratic was 
Beuudiiyea factor of the quartic. 

Wie Treserts of the variots@mrreducibility techniques 


are summarized as follows: 


SZ 





TABLE II. 


NUMBER WITH NUMBER OF 


LIMIT X DEGREE NUMBER OF NO LINEAR PReeevuCr BLE 
GRID MPWIC's FACTOR NPWIC's 
See 44 43 43 
Se aa JtgSi2 ZO 20 
> X 4 OH EON 96 


The results as illustrated in Table II complete the 
establishment of data bases of irreducible polynomials for 
three separate sets of polynomials upon which the utility 
of the developed transformational scheme and algorithms can 
meme omed. 9ihe irreducible polynomials for each grid are 


listed in Appendix A. 


53 





V. APPLICATION OF TRANSFORMATIONAL SCHEME TO DATA BASE 


imonmaenr  comacauire Knowledge of the effect ofthe 
transformational scheme on a given data base and hence ar- 
rive at some conclusions as to the utility of the scheme, 
the transformational scheme was applied to the data bases 
meseribed in Section IV, using a program called "'TESTER" 
which actually implements the testing of polynomial pairs 
Pmoecemibed in Seceron Ii]. The execution of TESTER ac- 
tually separates a set of polynomials into equivalence 
classes, each of which consists of all the polynomials of 
the original set which are related to each other through 
[fe transformation scheme. The application to the sets of 
polynomials resulted in the reduction of the 43 quadratic 
meneducible MPWIC's into 13 equivalence classes, 126 cubic 
iureducible MPWI€'s into 91 equivalence classes, and 96 
O@m@artic MPWIC's into 93 equivalence classes. These classes 
are tabulated in Appendix B. 

As a check on these results, an interactive program 
Gaited “XFORM’ was written and executed which actually 
performs desired sequences of D and T transformations upon 
meoclicchecspolynomual. Independently of the results of 
[TeawerR the same data bases were manipulated by the use 
of XFORM and the same results were achieved through trial 
and error transformation sequences. 

Also an interactive program called TESTERA was written 


moment 1s merely a minor revistonwor TESTER in that two 
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pelynomrals are entered via an interactive terminal and 
TESTERA determines whether or not a transformational rela- 


mronsiim exists and if so, describes the relationship. 


Se 





Vi) Rew, CONCLUSLONS, AND HYPOTHESES 


This research has shown that any polynomial which, when 
transformed to a monic polynomial, has integer coefficients 
and is CoMmemicit tOmtncne@mmdsrOrei1ts degree over the range 
of abcissae of the grid, can be tested for irreducibility 
by testing it for a transformational relationship to one 
mimper of Gach of the classes of irreducible classes of 
memeatcible MPWIG's for that particular degree grid. If 
a eneGisaShip exists with any of these classes then the 
feiynomial 1s irreducible. Note that the polynomial is 
meaced Tor anreducibility and is found to be reducible by 
‘default rather than the usual technique of attempting to 
determine reducibility and finding irreducibility by de- 
fault. | 

Since there are infinitely many polynomials which may 
Destested, the primary result of this research has been the 
establishment of an irreducibility criterion which can be 
easily and quickly pi het Eom ili idiate Seteot polyno- 
Mials. For instance, any quartic monic polynomial, say 
f(x), such that |£(x)|<5 for any 5 consecutive integer ar- 
guments can be tested for eeeedicdibadity by merely trans- 
forming the polynomial by means of the T transformation so 
that the S arguments are coincident to the abcissae of the 
degree 4, limit S grid. 

Deas beenotedt that, due to ge severely limited 


puslaswselected, the utility of the Detransformation is 
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demonstrable only for the data base for the quadratics. 

That is, the limit constraint of |f(x)|<5 precluded the 
mossibility of the existence of f£(x) and g(x) such that 

g(x) = qa” my aeand Moth etxyeand ¢(x) being coincident 
momenc orid tor depree greater than or equal to three. 
Mewever, in the quadratic case, this limit restriction does 
not preclude the possibility and results demonstrating the 
effectiveness of the D transformation were achieved. The 
execution of a modified version of TESTER which tested only 
ferme! transtormational relationships between members of the 
set of 43 irreducible quadratics resulted in the establish- 
Meme OL 16 Equivalence classes under this transformation. 

As previously noted, the testing for DT pairs of transforma- 
Ponseresulted in the establishment of only 13 equivalence 
classes. Comparison shows that 3 pairs of classes under 
mel transformation are compressed into 3 classes under the 
DiSthansrtormational pair. Thus, the utilization of the D 
transformation along with the T transformation extended the 
Se pabilitves OL the techniques. 

In view of the results of this work, pens ose arch 
icontie technique wollld seem warranted. This further re- 
search, if undertaken, could be to extend the technique in 
Sie sOmelcowe Of Several possible directions. Perhaps other 
and more powerful transformations could be identified and 
included imine transtLormatizonal seheme. Or modification 
OE Wewo Speeitications of the varios data sets might prove 
Sale Smet S auc noueas 1dentwaeds NO conceptual barrier 


Cepctwmer CxXtensSion. Merely relaxing the limit constraint 
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to some other fixed value, possibly the sum of the degree 
and some well-chosen constant, should provide interesting 
results when compared to the results of this work. Alter- 
natively, some variable limit such as f(x) = x" where n is 
the degree of £(X) might provide some additional insight 

into the total problem of a transformational scheme of ir- 


meauecibality criteria. 
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APPENDIX A 


PABRE 1 


NOMIALS 


ec if 


EGErPICTENTS LISTED AS 


AX2 + BX + € 
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AX4 + BX3 + CX2 + DX + E 


TABLE 3 


SMeETEICTENTS LISTED AS 
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| | tic | tieurdte I 1 ot I oe en | ee | itd 
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NOTES AND INSTRUCTIONS ON UTILIZATION OF PROGRAM "XFORM" 


io FORM" 2S ya pnecram whieh performs selected "T" and/or 
"D" transformations on a desired poiynomial. 

2. The degree of the polynomial is specified by the entry 
Sf aspositave integer. 

3. The coefficients of the desired polynomial are entered 
moeintegers separated by commas, representing the coeffi- 
ements of the polynommal an descendiamg powers of the vari- 
able. Zero coefficients must be included. 

4. The type of transformation to be performed is specified 
my entry of "Teor "D." 

be ihe tramslation factoueas entered as an integer. 

6. The dilation transformation is possible with an integer 
factor or a factor which is the reciprocal of an integer 
for certain polynomials. Thus, when the "D" transformation 
PmocrcGncd, anteger or recaprocal gtactoms must be specified 
imaecwcesponse of 1" or "RR." Then the actual factor is 
Seecrtted by the emtry of an integer value. If the factor 
as to be an integer, the value is the value entered. If 
i mitGeOm rs tO be a reciprocal the value is the récipro- 
cal of the integer value entered. 

7. The resulting polynomial is displayed as a series of 
integers separated by commas. These integers are the coef- 
frerencs of the descending powers of the variable. 

Seeeeene eal OPTION?” data entry pont allows the operator 


Mmsercer the Lesult Of the bast transtormation, the base 
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Ot the last transformation, the original polynomial for the 
present series of transformations, or the entry of a new 
polynomial as the next polynomial to be transformed by 
meats OL the appropriate response. The responses are "R," 
mee «6 0," and ""* respectively. 

9. Arithmetic data in the program is represented in binary 
integers of 31 digits. This allows a range of integers 
feom -262,144 to 262,143 for any single value. This re- 
striction is a limit of the IBM System/360 PL/1 implementa- 
LON . | 

10. The output coefficients are displayed to the operator 
as integers separated by commas. Internally, this is a 
Single character string constructed through various string- 
handiing operations. The maximum length of this string is 
Weecharactcrs including the commas, again due to the System 
/360 PL/1 implementations. 

11. The maximum number of coefficients for display as out- 
put is 64 due to problem data attributes in the program. 
Beever > Seino hestruction 15 insignificant in that the num- 


ber of coefficients is limited to 36 from Note #9 above. 
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load xform (xeq) 
EXECUTION BEGINS... 
DEGREE? 


= : 

COEFFICIENTS SEPARATED BY COMMAS 
mie Wale = 1,1 

TRANSFORMATION TYPE? 

= 

K = 


: 

io 7 6 LS, 19 

NEXT OPTION? 

_result ; 
TRANSFORMATION TYPE? 
ait 

K = 

a= 2 

1 el ee eo Mi ae | 

NEXT OPTION? 


SG 
TRANSFORMATION TYPE? 


ad 

A= "INTEGER" OR "RECIPROCAL"? 
_integer 

A = 

2 

1,°72,4,-8,-16, 32 

NEXT OPTION? 

Jf 

TRANSFORMATION TYPE? 

_d | 

A= "INTEGER" OR "“RECIPROCAL''? 
. 


_2 
eel, lye 1 
NEXT OPTION? 


> | 


DEGREE? 
“4 | 
COEFFICIENTS SEPARATED BY COMMAS 
ea een 
TRANSFORMATION TYPE? 
t 


~ | 
NO 


Mors 1495 51 
NEXT OPTION? 
(continued on next page)é 
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r 
TRANSFORMATION TYPE? 


_d 
A= “UNTEGER" OR “RECUPROCAL'? 


= 
a 


mee 
yulten 24, 992,596 
Neal OPTION? 


>| 


r 
TRANSFORMATION TYPE? 
d 
A= "INTEGER" OR "RECIPROCAL"? 
ak 
[RS 
a 
9, 517495,51 
NEXT OPTION? 
r 


RANSFORMATION TYPE? 
t 


COEFFICIENTS SEPARATED BY COMMAS 
i 0 Or | 
* TRANSFORMATION TYPE? 


ad 
A= "INTEGER" OR "RECIPROCAL"? 


>| 
WN 


BP) ye 
NEXT OPTION? 
O 


TRANSFORMATION TYPE? 
_d 
A= "INTEGER" OR "RECIPROCAL"? 


jf — 

we 

<1 005 —o4 
NEXT OPTION? 


DEGREE? 


END OF JOB 
Peer 6. 15 13.21.40 
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r 
TRANSFORMATION TYPE? 


_d 
A= "INTEGER" OR "RECIPROCAL"? 


i 


>| 


Bae 
lore 2,49 6 
NEXT OPTION? 


_r 
TRANSFORMATION TYPE? 


MUNVTEGER @OR “RECIPROCAL? 


Po 15 | 
NEXT OPTERON? 
r 


TRANSFORMATION TYPE? 


COEFFICIENTS SEPARATED BY COMMAS 
ieo0r 
* TRANSFORMATION TYPE? 


_d 
A= "INTEGER" OR "RECIPROCAL"? 
j 


>| 


O02 
XT OPTION? 


ze 
ms Ww 


lo 


ANSFORMATION TYPE? 


a 


"ENTEGER" OR "RECIPROCAL"? 


>|] 
—- 1 OW 


74 
> On Osos 
NEXT OPTION? 


>| 


DEGREE? 


END OF JOB 
Reem 2e1l/ 6.13 13, 2 bao 
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NOTES AND INSTRUCTIONS ON USE OF PROGRAM "TESTERA" 


1. “TESTERA" is a program which determines the transforma- 
tional relationship between two primitive polynomials with 
immereger coetfiveients. 

2. The "LIST" data entry specifies whether a single poly- 
nomial is to be tested against the first polynomial entered 
or a-succession of polynomials to be tested against the 
same polynomial. Response is 'Y' (yes)or 'N' (no). If 
Seis tT" is specified, Gre program returns to Ehe entry point 
for the second polynomial after completion of testing of’ 
each SOs nonilat 

3. The "DEGREE?" is entered as a positive integer. 

4. Coefficients cf polynomials are entered as integers 
separated by commas representing coefficients of descending 
soiers of variable. Zero coefficients must be included. 

S. Results of tests are displayed as integer values of K 
and A for transformations between the polynomials Or, 1G 

no transformational relationship exists, "NO TRANSFORMATION" 
1s Girsplayed. ~ 

Ceebrectcner @r both polynomials are Hon one ace they are 
transformed into monic polynomials prior to testing for 
melacionship between them, 

Pp oeoi2emand amemetientation restrretions for "“TESTERA” are 


the same as for "XFORM." 
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load testera (xeq) 
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_& 
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10, Gro 7 
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K=0...A=s3*%** 
DEGREE? 
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_yes 
DEGREE? 


2 
FIRST COEFFICIENTS SEPARATED BY COMMAS 
_1,1,1 
eae COEFFICIENTS 
ma DF) 
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NEXT C 
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Next mecOEFFICIENTS 


PEC ENTS 


m O D> 
iT] Of 
T1W 


ELST? 

_no 

DEGREE? 

cs 

PPRSot sCOEFFICERIDES SEPARATED BY COmMAS 
melee) la 1 1) 

New ecOmerlGt ENTS 

sig eptre, Lo, 52,04, 126525672, 10 Ae AUG 
K=0...A=2 

DEGREE? 


EesT? 


END OF JOB 
R: T=1.73/4.07 13.40.06 
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and discussed. 


The conclusion is drawn that a transformational scheme 
iceemiscrul tool for the investipationpote une mw a duet ane: 
of this type of polynomial and some hypotheses as to the 
Fetmmement and extension of the téehniquerarcems eaced. 
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